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^ Abstract 

pH , We extend a previous analysis of the lowest-lying octet baryon masses in covariant baryon chiral pertur- 

bation theory (ChPT) by explicitly taking into account the contribution of the virtual decuplet baryons. Up 
^ ■ to next-to-next-to-next-to-leading order (N'^LO), the effects of these heavier degrees of freedom are system- 



O 

m 



atically studied. Their effects on the light-quark mass dependence of the octet baryon masses are shown to 
be relatively small and can be absorbed by the available low-energy constants up to N^LO. Nevertheless, a 
better description of the finite- volume corrections of the lattice QCD data can be achieved, particulaiiy those 
with small M^L (< 4), which is demonstrated by a careful study of the NPLQCD and QCDSF-UKQCD 
small-volume data. Finally, we show that the predicted pion- and strangeness-baryon sigma terms are only 



X 

\ sUghtly changed by the inclusion of the virtual decuplet baryons. 
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I. INTRODUCTION 



In recent years, lattice chromodynamics (LQCD) simulations yi |2|] have made remarkable 
progress in studies of non-perturbative strong-interaction physics (e.g., see Ref. ^). Accurate 
computations of the lowest- lying octet baryon spectrum with nj = 2 + 1 dynamical simulations 
have been reported by several LQCD collaborations [|4|-jl2ll. Nonetheless, most LQCD calcula- 
tions employed larger than physical light-quark masses U and finite volume, therefore the obtained 
results have to be extrapolated to the physical point by performing the so-called "chiral extrapola- 



tion" 115 



1811 and taking into account the finite- volume corrections (FVCs) 11191 20 ] 



Chiral perturbative theory, as the low energy effective field theory of QCD [|2ll - l32ll . provides 



an appropriate framework to study the light-quark mass dependence and the lattice volume depen- 
dence of LQCD results. In the last decades, baryon chiral perturbation theory (BChPT) has been 
applied to study the ground state (g.s.) octet baryon masses with 11331 - 13611 and without H-B the 
explicit inclusion of the intermediate decuplet resonances. Recently, the = 2 + 1 LQCD simula- 
tions of the lowest-lying baryon masses at different combinations of light-quark and strange-quark 
masses have been studied with the SU(3) BChPT up to iHQ. It is found that the non- 

relativistic heavy baryon (HB) ChPT converges rather slowly up to next-to-next-to-leading order 



(NNLO) m 



4 in . The situation is much better in the finite range regulator (FRR) method P42!] and 



in the extended-on-mass-shel 



partial summation scheme [|45 



(EOMS) BChPT ill. Up to N^LO, the infrared BChPT [44], the 



47|] , and the EOMS BChPT MM all seem to be able to describe the 



LQCD data, but only in Ref. [|48l ] finite- volume effects were taken into account self-consistently 



and all the relevant LLCs were determined from the LQCD and experimental data. 
Using the covariant BChPT with the EOMS scheme 



studies on the LQCD mass results [143 . 



48 



SID. In Ref. 



150D, we have performed a series of 
431], it is shown that the EOMS BChPT 



can provide a better description of the PACS-CS Bg] and LHPCJ8D data and is more suitable for 
chiral extrapolation than the HBChPT up to NNLO. In Ref. Jsil], the NNLO EOMS BChPT is 
used to analyze the NPLQCD I12I] lattice volume dependent results at fixed quark masses. It 
is shown that for moderate Mt,L, FVCs need to be taken into account in order to describe the 



lattice data and they are also helpful to constrain the value of some relevant LLCs. In Ref. [|48 
we performed a simultaneous fit of all the publicly available LQCD data from the PACS-CS 



LHPC 1^, HSC M\, QCDSF-UKQCD Eil] and NPLQCD 

Up to now, only a few lattice collaborations have performed simulations at or close to the physical light-quark 
masses with fixed ~ mP'^^''- lll3lll4ll . 
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mi 



12 ] collaborations using the N^LO 



EOMS BChPT with the FVCs taken into account self-consistently. It is shown that the covariant 
SU(3) BChPT converges as expected with clear improvement order by order, and all the lattice 
simulations are consistent with each other, although their setups are quite different. It should be 
mentioned that the contributions of the virtual decuplet baryons were not explicitly included in 
the N^LO BChPT calculation of Ref. with the assumption that the decuplet effects can not be 
disentangled from those of the relevant LECs. 

However, in the SU(3) BChPT one should be careful about the contributions of the decuplet 
resonances since the average mass gap between the baryon octet and the baryon decuplet 6 = 
mo — rriQ ~ 0.3 GeV is similar to the pion mass and well bellow those of the kaon and eta 
mesons, Mk, Mrj. In Ref. [|33D, HBChPT was enlarged to include the decuplet and applied to 
calculate the octet baryon masses up to O(p^), and the importance of the spin-3/2 fields was 
pointed out. However, it was shown in Ref. [|34ll that the effects of the virtual decuplet on the 
octet baryon masses actually start out at 0{p^) in the same framework. For the spin-dependent 
quantities, the virtual decuplet contributions are found to be important in HBChPT, such as the 



SJ. In the EOMS BChPT, the effects 



magnetic moments [15211 and the axial vector form factors 
of the virtual decuplet are found to be negligible for the magnetic moments of the octet baryons 
if the "consistent" coupling scheme for the octet-decuplet-pseudoscalar coupling is adopted [15 511 . 
On the hand, up to NNLO, the virtual decuplet contributions seem to play an important role in 
descri bing the NPLQCD volume-dependent data [51] and in the determination of the baryon sigma 
terms [|56|1. Therefore, it is necessary to study the effects of the virtual decuplet on the light-quark 
mass and volume dependence of the LQCD data and on the determination of the baryon sigma 
terms at N^LO. 

In this work, we explicitly take into account the contributions of the virtual decuplet resonances 
to the g.s. octet baryon masses in the EOMS BChPT up to N^LO. The finite- volume corrections 
from the virtual decuplet are calculated self-consistently. Through a simultaneous fit of the pub- 
licly available ri/ = 2 + 1 LQCD data from the PACS-CS [i LHPC [8], HSC |9], QCDSF- 
UKQCD |1 11] and NPLQCD [12] collaborations, we perform a systematic study of the intermedi- 
ate decuplet effects on the chiral extrapolation of and the FVCs to the octet baryon masses, and on 
the determination of the octet baryon sigma terms through the Feynman-Hellmann theorem. 

This paper is organized as follows. In Sec. [Ill we collect the chiral effective Lagrangians 
involving the decuplet baryons and calculate their contributions to the octet baryon masses in the 
EOMS BChPT up to N^LO. In Seclml we perform a simultaneous fit of the LQCD data and study 



the effects of the virtual decuplet baryons in detail. A short summary is given in Sec. HVl 



n. THEORETICAL FRAMEWORK 



A. Chiral effective Lagrangians involving the decuplet baryons 

The baryon decuplet consists of a SU(3)-flavor multiplet of spin-3/2 resonances, which are 
represented with the Rarita-Schwinger field T"*^ = T^^^ (each element of T^'"^ is a four-component 
Dirac spinor). The physical fields are assigned to the tensor as T^^^ = T^^^ = A+/-\/3, 

T122 = A7v^, T222 = A-, Tl^3 ^ S*+/v^, = S*7v^, T223 = S^^/VS, = 

S*7v^, T233 = S*-/v^, and r333 = n-. 

The covariant free Lagrangian for the decuplet baryons is 

Ct = Tf' {t7'""'D^ - niDin ^ (1) 

where m/) is the decuplet-baryon mass in the chiral limit and D^T^^^ = d^T^^^ + (Fj,, T^Y^^, T^, 
being the chiral connection (see, e.g., Ref. ^\) and with the definition (X, T^)''^^ = (X);^T^^^ + 
{XY^T^^^+(X)^TI^'"^. In the last and following Lagrangians, we always apply the Einstein notation 
to sum over any repeated SU(3)-index denoted by latin characters a, b, c, ■ ■ ■ , and (X)^ denotes 
the element of row a and column b of the matrix representation of X. The totally antisymmetric 
gamma matrix products are defined as: 7^^^ = | [7^, 7^^], 7^^^" = |{7^'',7"} = —^£^'^'^"'^7/375, 
with the following conventions: gf^" = diag(l, —1, —1, —1), £0,1,2,3 = — £:°'^'2>3 = 1 and 75 = 
«7o7i7273- 

The C(p^) chiral Lagrangian for the decuplet baryons is: 

with x+ = 2% = 4i?odiag(m/, m^, nis) introducing the explicit chiral symmetry breaking, where 
mi and nis are the average light-quark and strange-quark masses. The parameters t^, tr, are two 
unknown LECs. 

Up to 0{p^) the chiral effective Lagrangian, describin g the interaction of the octet and decuplet 



baryons with the pseudoscalar mesons, can be written as [|55ll 

iC 



= -^e^'%d^Tf^)r'''B:d,<Pt + H.c, (3) 



(a) 



(b) 



(c) 



FIG. 1. Feynman diagrams contributing to the octet baryon masses with the intermediate decuplet reso- 
nances. The soUd Unes correspond to octet baiyons, the double lines to decuplet baryons, and the dashed 
lines denote pseudoscalar mesons. Black dots indicate an insertion from the dimension one chiral La- 
grangian (Eq. (|3]l), and black boxes (diamonds) indicate mass insertions. 



where we have used the so-called "consistent" coupling scheme for the octet-decuplet-pseudoscalar 



vertices 157 



58ll . The (p and B are the SU(3) matrix representations of the pseudoscalar mesons 



and of the octet baryons. The coefficient is the meson-decay constant in the chiral limit, and C 
denotes the (pBT coupling. 



The propagator of the spin- 3/2 fields in d dimensions has the following form ll59ll 



f + mr, 



m|) + ie 



1 



1 



{d — l)m£i 



d-2 



{d — l)m 



D 



(4) 



B. Virtual decuplet contributions to the octet baryon masses 

Because the baryon mass, which is of the same order as the chiral symmetry breaking scale 
AchPTj does not vanish in the chiral limit, a systematic power-counting (PC) is destroyed 
beyond the leading order calculation in BChPT [|24|l. In order to restore the chiral power- 
counting, the extended-on-mass-shell (EOMS) renormalization scheme was proposed [49, Isoj . 
The essence of the EOMS scheme is to perform an additional subtraction of power-counting 
breaking (PC B) p ieces beyond the MS or MS renormalization scheme. Different from the infrared 
(IR) BChPT ieOD and HBChPT |I6l|], the EOMS BChPT is not only covariant, but also satisfies 



all ana^ 



faster pS 



yticity 



63 



and symmetry constraints (see, e.g., Ref. [|62|] ). In addition, it converges relatively 



64|l . In this work we use the EOMS scheme to remove the PCB terms from the 



one-loop diagrams. 

The octet baryon masses up to N^LO and with the virtual decuplet contributions can be written 

as 



niB 



, (2) , (3) , (4) 

mo + + + 



m 



(D) 
B ■ 



(5) 



TABLE I. Coefficients of the NNLO virtual decuplet contribution to the self-energy of the octet baryons 
(Eq. ®). 



A S H 

16/^2 /ui2 8^2 4/^2 

?_B,7r 3 ^ 9*" 3*^ 

Ja) 4^2 8/^2 40/^2 /f/^2 
3*^ 3*^ 9 ^ 

Ao) n n 4/^2 4/^2 



TABLE II. Coefficients of loop diagrams (Fig.[TIb/c)) to the self-energy of the octet baryons (Eq. ([6)). 



N 
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Ab/c) _ 




p(b/c) 
?AS*, 


Ab/c) 


= c' 


Ab/c) 



q/^2 _ 8/^2 _ ^2 

"^•^ ''SA,/^ ~ 3*^ ^SS*,K — ^ 

- 9(^2 _ 2/^2 tCVc) _ r>2 

Ab/c) _ p2 Ab/c) _ p2 

Ab/c) _ 2p2 Ab/c) _ „p2 

sss*,x — 3'" ^e:a-,x ~ 



TABLE III. Coefficients of LO contribution to the self-energy of the decuplet baryons (Eq. (H)). 

A S* H* VL- 



£,D,TT to + 3tD tQ + t£) to — to to — 3tD 



Here, mo is the chiral limit octet-baryon mass, and m^^\ m^^\ m^"* correspond to the 0{p'^), 
0{p^), 0{p^) contributions from the octet-only EOMS BChPT, respectively. Their explicit ex- 



pressions can be found in Ref. ll48ll . where the PCB terms are already removed. The last term 



m^j^^ denotes the contributions of the virtual decuplet resonances up to N'^LO. After calculating 
the Feynman diagrams shown in Fig. [Hand subtracting the PCB terms with the EOMS scheme, 
the virtual decuplet contributions to the octet baryon masses can be expressed as 



D=A, s'*, H*, n- 

The first term of Eq. (|6]) is the NNLO contributions of Feynman diagram Fig. [IJa). The loop 



function Hf{M^) can be found in Ref. [|43ll and the corresponding coefficients are listed 
in Table m The next term is the virtual decuplet contribution at 0{p^) from the one-loop diagram 
of Fig. [TJb) and the wave function renormalization diagrams of Fig. [IJc). The Clebsch-Gordan 
coefficients .^bd!/, are tabulated in Table HIl and the loop function H^^^^\m^) has the following 
form: 

1 



where W 



24mQm|| 



m 



(2) 



Ml [2(1 Im^ + SmoTUD + 9m^) 



x{ml - ml) + (5m2 + Smom^ + 18m^)M| - 6Mj] 



mg^M| [-34m^ - 24m|jmB + 30m^ 



SGniQml 

-6m^Mj - 6Mj + 6momD{Aml + Mj) + 3m2(4m^ + 7Mj)] 



— ^mg^Mjln^ [l2ml,(mo + mz,) 

+6(mo - rnlml - 2mom% - 3m%)Ml + 4(mo + rnQmo + 3m|,)M^ 

1 O^. /A'/^X 



6mQm^ 



^mg)M| In 



[m|,(9mo + 8m£)) 



+3(mg - ml){2ml + moms + 2m^)M| - mo(4mo + mB)Mj + 



12mQm|5 



M^(mo - mo) (mo + m^) 



In 



m/) (— 5mQ + 2mQmD — 3m^) + 2mo (mg — m^mr) + 3m|,) m 



2 2 

mj, — mQ 



,(2)' 
D 



6m^ 



moMMn 



V ) 



('2) (2) 

6m£)(mo + ra£))m^ — mo(4mo + 'ini£,)m\) 



rrin 



2mQm£) + m^, — M^)(mQ + 2mQmD + — 



2\2 



(2) , 2 (2) 
D +^D^B 



,{2) 



,(2) 



—niQmD i^rn^B + 2m~'^\ + 2mQm^'' 



arctan 



m 



D 



- Mi 



+ arctan 



m 



D 



Ml 



(7) 



—m, 







I) — M^)^ + 2ml{m1 + M^), the NLO octet baryon mass m^^ is given 



in Eq. ^ and the NLO decuplet baryon mass is 



m 



(2) 
D 



(8) 



The corresponding coefficients axe listed in Table Hill It should be noted that in obtaining 



the results of Eq. the decuplet-octet mass difference, 5 = nio — mo, is considered up to all 
orders n43n . 



C. Finite-volume corrections from the virtual deucplet baryons 



To study the 



LQCD data from the PACS-CS \M, LHPC m, HSC 



10], 



QCDSF-UKQCD lUO 



and NPLQCD [|l2h collaborations, FVCs should be taken into account. That's because at present 
all LQCD simulations are performed in a finite hypercube with lattice sizes L ~ 3 — 5 fm and 
the finite-volume effects on the simulation results cannot be neglected. Recently, the FVCs to the 
octet baryon masses have been studied in the HBChPT ll65n and the EOMS BChPT 115 ill up to 
NNLO. It was shown that the FVCs can be helpful to constrain some relevant LECs. In order to 
self-consistently take into account the FVCs to the octet baryon masses, we need to calculate the 
FVCs up to N^LO in the EOMS BChPT. Following the same procedure as detailed in Ref. 115 ill , 
one can easily calculate the C(p^) BChPT results in a finite hypercube by replacing the H of 
Eq. ([5]) with H = H + 6G. For the purpose of this work, we only present the FVCs from the 



virtual decuplet baryons. The octet-only formulas of FVCs can be found in Ref. [|48 
The FVCs to the loop results of Fig. [TJa,b,c) are 



6G 



and 



(a) 



dx 



ml (mo(l - x) + nir)) 



{Ml) 



ml (mo(l — x) + mo) M. 
6ml 



^^3/2 (Ml) 



(9) 



6G 



{b/c)_ mp 



B,D 



dx 



+ 



2ml{x — l)m 



3mlmDx{x 



D^m^"* — 2m^AiKx — l)m^^' + 3mQm^x(x 



jj — m^mr) ( + Sm^-* {x — 1] 



.(2) 



X 6i/2 {Ml) 



m 



(2) 
D 



m 



+momDMl ( 3m^^\x — 1 



m 



{'2') ■ ■ z 

- "im^m'^^ml X — 2m^A^^m^ 
momo^x 



(2) 



+3momDMlx —ml{x 



X ^3/2 {Ml) 



T \ I (2) (2) 



2 (2) 



X55/2 {Ml)} 



(10) 



respectively. Here Ml = x'^ml — x{ml — ml) + (1 — x)M^ — ie, and the "master" formulas are 
defined as 



r,-l/2-r( fT72\3-2r 



TT 



3/2r(r) 



where i^^ri (2;) is the modified Bessel function of the second kind, and ^ = Yl Yl (1" 
6{\n\,0)) withn = {n^,ny,n-^). 



-oo n,=— oo 



ni. RESULTS AND DISCUSSIONS 



In this section, the effects of the virtual decuplet baryons on the g.s. octet baryon masses are 



systematically studied by fitting the Uf 

n " 



2 + 1 LQCD data of the PACS-CS M, LHPC 
HSC m, QCDSF-UKQCD HH] andNPLQCD |Q] collaborations. 



A. Light-quark mass dependence of the octet baryon masses 



Up to N^LO, there are 19 unknown LECs (mo, h^, hjj, hp, hi^^, and di_^^ 7 §) in the octet- 
only EOMS BChPT. To take into account the contributions of the decuplet baryons, one has to 
introduce four more LECs, m^), t^, and C (see Section IIA). The (pBD coupling constant C can 
be fixed to the SU(3)-average value among the different decuplet-to-octet pionic decay channels, 
i.e., C = 0.85 ll56n rl A moderate variation of C has no significant effects on our final results. The 



LECs to, tn, and m/5 can be fixed by fitting the NLO decuplet mass formula Md = mo — rri^o 
the physical decuplet baryon masses. Because to and rriD cannot be disentangled at the physical 
point, one only obtains a combination of mo and to with m^f = ttid — to{2Mj^ + M^) = 1.215 
GeV and to = —0.326 GeV~^ In the following, the octet-decuplet mass splitting 5 = mo — rrio 
is fixed to be 0.231 GeV-the average mass gap between the octet and decuplet baryons. Therefore, 
one can fix the four LECs in the following way: mo = mo + 0.231 GeV, to = (mo — 0.984) / 0.507 
GeV"\ to = -0.326 GeV"^ and C = 0.85. As a result, the same 19 LECs as those in the octet- 



only BChPT need to be determined. The other coupling constants are fixed as in Ref. [14811 : the 



meson decay constant = 0.0871 GeV, the baryon axial coupling constants D = 0.8, F = 0.46 
and the renormalization scale = 1 GeV. 



Following the same procedure as in Ref. [|48ll . we use the formulas Eg. jSl) to fit the octet 
baryon masses of the LQCD simulations and of their experimental values 116611 . The LQCD data 
to be studied are taken from the PACS-CS M, LHPC M, HSC m, QCDSF-UKQCD lul] and 
NPLQCDJIJL data satisfying AU < 500 MeV and M^L > 4, named as Set-I in Ref. MM- 
5511 the value of C is fixed from tiie A(1232) nN decay rate, which yields C = 1.0. But in our 



In Refs. 143 



previous study of the NPLQCD data, this couphng turned out to be somewhat smaller pill 



9 



The restrictions on the lattice data are taken to ensure that the N^LO BChPT is valid for these 
pion (light-quark) masses and lattice volumes. Later we will slightly relax the restrictions to test 
whether the inclusion of the virtual decuplet baryons can extend the applicability region of the 
BChPT. The so-obtained values of the LECs from the best fit and the corresponding x^/d.o.f. are 



tabulated in 



able|IVl For comparison, we also list the octet-only best fit results of lattice data Set- 



I from Ref. HaA . From the x^/d.o.f., one can conclude that the inclusion of the virtual decuplet 
baryons does not change the description of the LQCD data. On the other hand, the values of the 
LECs have changed a lot, as can be clearly seen from Table lIvR This confirms the assumption that 



using only the octet baryon mass data, one can not disentangle the virtual decuplet contributions 



from those of the virtual octet baryons and the tree-level diagrams [|48l1 . In other words, for the 
static properties of the octet baryons, most contributions of the virtual decuplet are hidden in the 
relevant LECs, as one naively expects. Below, we will see that their inclusion, however, does 



improve the description of the volume-dependence of the LQCD data, as also noted in Ref. 15 111 . 

In Fig. |2l setting the strange-quark mass to its physical value, we show the pion mass depen- 
dence of the octet baryon masses in the N^LO EOMS BChPT with and without the virtual decuplet 
baryon contributions. It is clear that the two N^LO fits give the same description of lattice data 
Set-I, as can be inferred from the same x^/d.o.f. shown in Table HVl 

B. Finite-volume corrections to the octet baryon masses 



In Ref. US in . we have studied the FVCs to the g.s. octet baryon masses using the EOMS BChPT 
up to NNLO, and found that the finite-volume effects are very important and cannot be neglected. 
Therefore, in this work the FVCs are self-consistently included in Eq. © to analyze the lattice 



data. The NPLQCD [|12|] simulation is performed with the same pion mass of Mj, ^ 390 MeV 
and at four different lattice sizes L ~ 2.0, 2.5, 3.0 and 3.9 fm. Therefore, it provides a good 
opportunity to study the FVCs to the octet-baryon masses. 

In Fig. [H we contrast the NPLOCD data with the N^LO EOMS BChPT using the Fit-I LECs 



from Table lYl As stated in Ref. MM, three sets of the NPLQCD data with M^L > 4 are included 
in lattice data Set-I and denoted by solid points in Fig.[3l Another set with M^L = 3.86 (hollow 
points) is not included. Both the octet-only and the octet plus decuplet (0-i-D) BChPT can give 



^ The same phenomenon has been observed in the studies of the octet baryon magnetic moments II55I1 and the octet 
baryon masses up to NNLO ll56ll . 
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FIG. 2. (Color online). Pion mass dependence of the LQCD data in comparison with the best fits of the 
EOMS BCiiPT up to N^LO with (solid lines) and without (dashed lines) the virtual decuplet contributions. 
The lattice data have been extrapolated to the physical strange-quark mass and infinite space-time. 

a reasonable description of the FVCs. At the e~'"''^/(m^L) < 0.2 region, these two fits give 
essentially the same results. With the increase of e~'"''''^ / {ni-^L) ( the decrease of lattice size L), 
the 0-i-D BChPT results are in better agreement the NPLQCD data, especially for the nucleon 
mass. It seems that the virtual decuplet baryons can help to improve the description of the FVCs, 
although the BChPT results are still a bit larger than the LQCD data at small M^L. 

It is interesting to check whether the 0-i-D best fit can describe the lattice data with larger pion 
masses and/or smaller lattice volumes. In Fig.gl the PACS-CS, LHPC, HSC and QCDSF-UKQCD 
lattice data with < 700 MeV are compared with the best N^LO 0-i-D EOMS BChPT with the 
Fit-I LECs of Table HVl The lattice points included in the fit are denoted by solid points and those 
excluded in the fit by hollow points. It is clear that the N^LO BChPT can describe reasonably well 

11 



the LQCD data, even those excluded in the fit. The average deviation of the BChPT results from 
the LQCD data, defined asB 

1 ^LQCD / j^i _ j^i \ 2 

~2 _ I ^'^LQCD ^^^BChPT 

A^LQCD Y ^LQCD 

is 3.1, 2.5, 1.2 and 1.2 for the PACS-CS, LHPC, HSC, and QCDSF-UKQCD data, respectively. 
Here, it should be noted that in Fig. H only the QCDSF-UKQCD data with A^^ = 32 are shown 
and those simulated in a smaller volume with Ng = 24 are not explicitly displayed. Including 
them in the x^, one would have obtained a = 22.3. 

It is clear from the above comparisons that using the LECs determined from the best fit to 
lattice data Set-I, the BChPT cannot well describe the LQCD data obtained in smaller volumes, 




1.14 I — ' — I — ' — I — ' — I — ' — I — ' — I — ' — I 
0.001 0.002 0.003 0.004 0.005 0.006 

Exp[-mjj L]/m^L 




Exp[-mjj L]/mjjL 




FIG. 3. (Color online). Lattice volume dependence of the NPLQCD data in comparison with the EOMS 
BChPT up to N'^LO with (solid lines) and without (dashed lines) the virtual decuplet contributions. The 
three black points with M^L > 4 are included in data Set-I, while the hollow points with M^L = 3.86 are 
not. 



The uncertainty of the lattice data, A^qcd, can be found in Ref. 14811 . 
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TABLE IV. Values of the LECs from the best fit to the LQCD data and the experimental data at 0{p'^) (see 
text for details). 





Fit-I 


Fit 1 9 


Fit 1 ^ 




w/o decuplet 


w/ decuplet 


w/ decuplet 


w/ decuplet 


iiiQ [ivie V J 


879(22) 


908(24) 


Q1 ni' 1 7"! 


Qi n('9'^'i 




-0.609(19) 


-0.744(16) 


— u. ( 1 ^10 j 


n 779M R^i 


u D [\jc V J 


0.225(34) 


0.355(20) 


n '^^9('9i 


n '^^'^/'9n'\ 

U.OfJtJ \^ZU j 


U p [VJC V J 


-0.404(27) 


-0.552(28) 




— U.OOOt^ZD J 


0\ \\jc V J 


0.550(44) 


1.08(6) 


L.oy\ 1 ) 


1 ^9('S'l 
i.OZ(^o J 


02 L*-J£ ^ J 


-0.706(99) 


0.431(93) 


U.DOU\^ ( Z j 


U. 000(^100 j 


03 [vje V J 


-0.674(115) 


-1.83(15) 




1 Qf^n 7^ 


64 [GeV~^] 


-0.843(81) 


-1.57(4) 


-1.67(2) 


-1.59(6) 


t/5 [VJC V J 


-0.555(144) 


-0.355(74) 




1 'i9('?\ 
— 1 .OZl^Z J 


[VJC V J 


0.160(95) 


-0.423(117) 




— u.zy / t^oo J 


1/7 [VJC V J 


1.98(18) 


2.79(15) 




9 fi'^/'i^'i 

Z.DO\^0 J 


[vjc V J 


0.473(65) 


-1.73(6) 




1 Qfi^'9'\ 
— i.yui^z ^ 


ui [oe V J 


0.0340(143) 


0.0157(130) 




n nn/i 1 ^^^^ 


U2 [oe V J 


0.296(53) 


0.445(57) 




U. 404(^00 j 


da [GeV~^] 


0.0431(304) 


0.328(18) 


0.280(22) 


0.264(15) 


(7^ rr;pv~3n 


0.234(67) 


-0.117(59) 


—0 0345('676l 


n no5Qni'5393'i 

\J .\J\J0cf\J\00c/0 J 


d5 [GeV-3] 


-0.328(60) 


-0.853(77) 


-0.831(59) 


-0.883(69) 


d7 [GeV-3] 


-0.0358(269) 


-0.425(39) 


-0.464(12) 


-0.497(33) 


ds [GeV-3] 


-0.107(32) 


-0.557(56) 


-0.602(18) 


-0.651(50) 


xVd.o.f. 


1.0 


1.0 


1.0 


1.2 
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FIG. 4. (Color online). The PACS-CS, LHPC, QCDSF-UKQCD and HSC lattice data in comparison 
with the 0+D BChPT best fit as functions of the pion mass. The lines in each panel (from bottom to 
top) correspond to N, A, S and H, respectively. The kaon mass is fixed using Mj^ = a + bM^ for the 
con^esponding lattice ensemble with a and b determined in Ref. [48;]. The lattice data have been extrapolated 
to infinite space-time using the corresponding BChPT fit. X-,^ = ^ (M^ + 2M^)/3, Xn = {niN + m-^ + 
ms)/3, where the meson and baryon masses are the physical ones. 



particularly those of the QCSDSF-UKQCD data with Ns = 24. On the other hand, the virtual 
decuplet contributions seem to be helpful in this regard. Furthermore, it should be noted that we 
have chosen lattice data set-I by requiring < 500 MeV and M^^L > 4. These criteria yielded 
a x^/d.o.f. = 1, but nevertheless, are a bit arbitrary. In the following subsection, we would like 
to slightly relax the above criteria and study whether the LQCD data with smaller Mt^L can be 
described at a reasonable sacrifice of the x^/d-o.f .. 
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FIG. 5. (Color online). Lattice volume dependence of the NPLQCD data in comparison with the best fits of 
the N'^LO 0+D BChPT. The solid lines correspond to the results of Fit-I, the blue-dot-dashed lines to the 
results of Fit- 12 and the green-hashed lines denote the results of Fit- 13. 



C. Description of the NPLQCD and QCDSF-UKQCD small-volume data 

In order to better describe the LQCD data obtained in smaller volumes, we have to take into ac- 
count in the fit a few more LQCD data with small M^L. Furthermore, to guarantee that the N'^LO 
BChPT is still valid, we have to ensure that the corresponding x^/d.o.f. be around 1. Therefore, 
we slightly relax the criterium of M^L > 4 to M^L > 3.8 and keep that of < 500 MeV. 
Fourteen sets of LQCD data from the PACS-CS, LHPC, HSC, QCDSF-UKQCD and NPLQCD 



collaborations satisfy the new criteria Hm . 



Because the NPLQCD results 111211 are obtained with the same light-quark masses but at differ- 
ent lattice sizes, they are more suitable for studies of FVCs and can better constrain the relevant 
LECs. First, we add the NPLQCD data point with M<^L = 3.86 to lattice data Set-I to form a new 
data set-Set- 12, and fit it using the N^LO 0-i-D BChPT. In Table lYl the LECs from the best fit 

15 




FIG. 6. The QCDSF-UKQCD lattice data in comparison with the N^LO 0+D BChPT results of Fit- 12 (left 
panel) and Fit-13 (right panel). The black points are included in the data Set-12/13, the red (with Ng = 24) 
and blue (with Ns = 32) points are not. The FVCs have been subtracted from the lattice data using the 
corresponding BChPT fit. 

(named Fit-12) and the corresponding x^/d.o.f. = 1.0 are listedlfi Compared with the 0+D Fit-I 
results, most of the LECs do not change much, particularly those of bo, ho, hp and mo- In Fig. [51 
the finite- volume dependences of the NPLQCD data are compared with the N^LO 0-i-D BChPT 
results calculated using the Fit-12 LECs. The BChPT results are in very good agreement with 
the lattice data and the corresponding is 0.7. It seems that the Fit-12 LECs can give a better 
description of the LQCD data with small M^L. This is further confirmed by contrasting the Fit-12 
results to the QCDSF-UKQCD lattice data. As shown in Fig.[6l both the data with Ns = 32 (green 
points) and Ng = 24 (red points) are better described than the BChPT results using the Fit-I LECs. 
The corresponding are 0.6 and 7.3, respectively. 

We note that the NPLQCD simulation is performed at fixed M^^ ^ 390 MeV and close to the 
physical strange-quark mass, and the values of M^L, M^L are all larger than 5 (for the smallest 
lattice size L = 2.0 fm, MkL = 5.4 and M„ L = 5.8). On the other hand, the situation is different 



for the QCDSF-UKQCD collaboration [|l IQ . They start with a common sea quark mass near 
{2m^^' + m^^^^-)/3 and split the masses symmetrically n67n . Thus, for the data with Ng = 24, 
not only A/^L is smaller than 4, but MkL and M^L are all around 4. 

In order to take into account the finite- volume effects induced by the smaller MkL and Mr^L, 
we add to data Set-12 the LQCD data set with 3.8 < M^L < 4 from the QCDSF-UKQCD 



If the octet-only BChPT were used to fit data Set-12, the best x^/d.o.f . would be 1.1. 
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collaboration. The new data set-Set- 13-has 13 lattice data points. Using the N^LO 0+D BChPT 
to fit Set-13, we obtain the results tabulated in Table HVl (named as Fit-13). Compared to Fit-12, the 
values of the LECs do not change dramatically, and the corresponding x^/d.o.f. becomes 1.2u- In 
Fig. ID and Fig. |6l we show the N^LO 0+D BChPT resuks calculated using the Fit-13 LECs. It is 
clear that the description of the NPLQCD data becomes a bit worse (x^ = 1.1), while the QCDSF- 
UKQCD data can be better described, with = 0.6 for A^^ = 32 and = 4.2 for A^^ = 24 lattice 
data, respectively. Nevertheless, it seems that up to N'^LO, the 0-i-D EOMS BChPT still cannot 
achieve a perfect description of the QCDSF-UKQCD data with A'^, = 24. 

It should be pointed out that because Fit-12 can describe better the volume dependence of the 
NPLQCD and QCDSF-UKQCD data while still yields a xVd-o.f. = 1.0, it is preferable than 
Fit-13 and Fit-I. 



D. Pion- and strangeness-baryon sigma terms 



It was pointed out in Ref. [|56D that the virtual decuplet baryons may increase the predicted 
pion-nucleon and decrease the strangeness-nucleon sigma terms. In addition, in Ref. [|48|] . the 
predicted a-^N is indeed smaller than that of the NNLO 0-i-D BChPT [ | 43i1 while the asN is larger. 
It is interesting to check whether similar effects still exist at N^LO. 



TABLE V. Sigma terms of the octet baryons predicted by the N'^LO BChPT with LECs determined from the 
global fits to different lattice data sets. The first error is statistical and the second is systematic, estimated 
by taking half the difference between the N'^LO result and the NNLO result. 







O'sB 


Fit-I 


Fit-12 
w/decuplet 


Fit-13 
w/decuplet 


Fit-I 


Fit-12 
w/decuplet 


Fit-13 
w/decuplet 


w/o decuplet w/ decuplet 


w/o decuplet w/ decuplet 


iV 


43(1)(6) 46(2)(12) 


47(1)(12) 


47(2)(12) 


126(24)(54) 157(25)(68) 


149(22)(63) 


157(22)(61) 


A 


19(1)(7) 20(2)(13) 


21(2)(13) 


22(2)(12) 


269(23)(66) 256(22)(60) 


250(25)(54) 


257(19)(51) 




18(2)(6) 19(2)(6) 


20(2)(7) 


21(2)(6) 


296(21)(50) 270(22)(47) 


266(23) (46) 


272(20)(50) 




4(2)(3) 6(2)(5) 


7(2)(4) 


7(2)(5) 


397(22)(56) 369(23)(50) 


366(23) (48) 


372(21)(49) 



If the octet-only BChPT were used, the obtained x^/d.o.f . would be 1.7. 
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In Table |Vl we tabulate the predicted a^B and cr^s by the N^LO 0+D BChPT with the Fit- 
I, Fit- 12, and Fit- 13 LECs using the Feynman-Hellmann theorem. The only-octet results from 
Ref. [i48i1 are also listed for the sake of comparison. The inclusion of the virtual decuplet baryons 
only slightly increases the pion-baryon sigma terms and the strangeness-nucleon sigma term and 
decreases the other strangeness-baryon sigma terms. 

We have checked if a smaller strangeness-nucleon sigma term, such as 21 ±6 MeV as predicted 



in Ref. 116 811 . were used as a constraint, the x^/d.o.f from a simultaneous fit of data Set-I and the 



small asN would increase to about 3. In Ref. ll69ll . using the R^.CS-CS data as an example, the 



influence of scale setting in LQCD simulations on the values of sigma terms is discussed. A 
comprehensive study of all the other LQCD data may be needed to understand the effects of scale 
setting on the results of a global study such as that performed in the present work. 



IV. SUMMARY 



We have calculated the virtual decuplet contributions to the ground state octet baryon masses 
within the EOMS BChPT up to N'^LO. Finite-volume corrections are calculated self-consistently. 
Through a simultaneous fit of the latest nj = 2 + 1 LQCD simulations from the R\CS-CS, LHPC, 
HSC, QCDSF-UKQCD and NPLQCD coUaborations, the effects of virtual decuplet baryons on 
the light-quark mass and volume dependence of the LQCD data are systematically studied. 

It is shown that the contributions of virtual decuplet baryons affect little the light-quark mass 
dependence of the octet baryon masses, indicating that their effects can not be easily disentangled 
from those of the virtual octet baryons and the tree-level diagrams. On the other hand, a slightly 
better description of FVCs can be achieved once the virtual decuplet baryons are taken into ac- 
count, especially for the lattice data with M^L < 4. This is demonstrated by a careful study of the 
NPLQCD and QCDSF-UKQCD small-volume data. 

Regarding the pion- and strangeness-baryon sigma terms of the octet baryons, it is shown that 
at N^LO, the effects of the virtual decuplet baryons are small. Effects of lattice scale setting may 
need to be studied to understand the relatively large asN predicted in the present work. 
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